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Bulge growth through disk instabilities in 
high-redshift galaxies 

Frederic Bournaud* 


Abstract The role of disk instabilities, such as bars and spiral arms, and the asso¬ 
ciated resonances, in growing bulges in the inner regions of disk galaxies have long 
been studied in the low-redshift nearby Universe. There it has long been probed 
observationally, in particular through peanut-shaped bulges (Chapter 6.3). This sec¬ 
ular growth of bulges in modern disk galaxies is driven by weak, non-axisymmetric 
instabilities; it mostly produces pseudo-bulges at slow rates and with long star- 
formation timescales. Disk instabilities at high redshift (z >1) in moderate-mass 
to massive galaxies (lO*** to a few lO" M© of stars) are very different from those 
found in modern spiral galaxies. High-redshift disks are globally unstable and frag¬ 
ment into giant clumps containing 10^^^ Mq of gas and stars each, which results in 
highly irregular galaxy morphologies. The clumps and other features associated to 
the violent instability drive disk evolution and bulge growth through various mech¬ 
anisms, on short timescales. The giant clumps can migrate inward and coalesce 
into the bulge in a few 10® yr. The instability in the very turbulent media drives in¬ 
tense gas inflows toward the bulge and nuclear region. Thick disks and supermassive 
black holes can grow concurrently as a result of the violent instability. This chap¬ 
ter reviews the properties of high-redshift disk instabilities, the evolution of giant 
clumps and other features associated to the instability, and the resulting growth of 
bulges and associated sub-galactic components. 


1 Introduction 


High-redshift star-forming galaxies mostly form stars steadily over lo ng timescales, 
merg er-driven starbursts being only a minority of galaxies (chapter 6.1 : lLe Floch et al 


rgei 

2005h . At redshifts z > 1, moderate-mass and massive star-forming galaxies (10 


10 


to a few 10** Mq of stars) have rapid gas consumption timescales and stellar mass 
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doubling timescales, of the order of a Gyr at z = 2, depending mostly on red- 
shift and weakly dep ending on mass, with rare deviations to the mean timescale 
( Schreiber et all2014 ). 

These star-forming galaxies have very irregular morphologies in the optical, es¬ 
pecially compared to nearby spirals disks of similar mass, as unveiled by deep sur¬ 
veys over the last two decades. They also have very high gas fractions, about 50% of 
their baryonic mass, as probed recently with interferometric studies. The high gas 
fractions and mass densities cause strong gravitational instabilities in the galactic 
disks, which results in disk fragmentation, and causes very irregular, clumpy mor¬ 
phologies. These irregular morphologies are often dominated by a few giant clumps 
of 10^^® Mq of baryons, rotating along with the host galaxy. 

This violent instability can drive the formation and growth of bulges, either by 
inward migration and central coalescence of the giant clumps and/or by gravitational 
torquing of gas and instability-driven inflows. This Chapter reviews the properties of 
the violent instabilities, and the pieces of evidence that the clumpy morphologies are 
caused by such violent disk instability rather than mergers or other processes. It then 
reviews the evolution of the giant clumps, their response to intense star formation 
and associated feedback processes, and the properties of bulges formed through this 
process. It eventually reviews recent results on other sub-galactic structures (such 
as thick disks and central black holes), which may grow concomitantly to bulges 
through high-redshift disk instabilities, and compares the role of this high-redshift 
violent instability to the contribution of low-redshift secular evolution through weak 
instabilities such bars and spiral arms. 


2 Clumpy galaxies and the violent disk instability at high redshift 


This section reviews the properties of star-forming galaxies at high redshift, the 
observed signatures of the underlying disk instabilities, and the main related theo¬ 
ries. Throughout, we consider galaxies at redshift z ~ 1—3, with stellar masses of 
lO'** to a f ew 10*^ Mf7i, and that are “normally” s t ar forming on the so-c alled Main 
Sequence ( Daddi et al. . 20071 Elbaz et al. . 2011 : Schreiber et al. . 2014I) . i.e. with 
specific star formation rates of the order of a Gyr^* at redshift z = 2, as opposed to 
rare starbursts with faster star formation. 


2.1 Clumpy galaxies at redshift 1-3: global morphology 


The characterisation of the structure of star-forming galaxies at high redshift has 
steadily developed over the last two decades, driven mainl y by deep surveys from 


the Hubble Space Telescope in the optical wavelengths (e.g. lCowie et al.Lll996l) . and 


later-on in the near-infrared, but also accompanied by modern techniques to iden - 
tify and select star-forming galaxies in these deep surveys (e.g. Daddi et al.Ll2004h . 
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The highly irregula r structure of high-redshift galaxies was first pointed out in the 


Hubb le Deep Field dAbraham et al.l 1 19961: Ivan den Bergh et aU 1 19961: ICowie et^ 


19961) . Star-forming galaxies appeared to have highly irregular morphologies, dom¬ 
inated by a few bright patches, with the striking example of the so-called “chain 
galaxies” where the patches are almost linearly aligned. While reminiscent of 
nearby dwarf irregulars, these morphologies where found in galaxies 10-100 times 
more massive - a mass regime at which nearby galaxies are almost exclusivel 


regular disk-domi nated galaxies, most often barred spirals (lEskridge et al.L 1200' 


i 


Block et al.L 120021) . or spheroid-dominated early type galaxie s. The lack of regular 


barred spirals, suspected in s uch deep optical imaging surveys (Ivan den Bergh et al 


19961: Abraham et al. . 1999), actually required deep-enough near-infrared surveys 


to be confirmed: otherwise the irregular structure could result from band-shifting 
effects, namely the fact that optical observations of z > 1 objects probe the ultra¬ 
violet emission, strongly dominated by young star-forming regions, rather than the 
underlying mass distribution dominated by older stars. The first near-infrared sur¬ 
veys unveiled counterexamples of high-redshift gal axies with a more regular disk 


structure in the underlying older stellar populations (ISheth et al.L 120031) . yet deeper 


and wider fields eventually c onfirmed the gradual disappearance of regular (barred) 


spiral disks atz > 0.7 — 1.0 ( Sheth et"^ 2008t Cameron et al. . 2010l : 
2014; Simmons et al.L 2014 ). 


Melvin et al 


Even before infrared data could resolve kiloparsec-scale structures at z = 2, de¬ 
tailed spatially-resolved stellar population studies were used to reconstruct the stel¬ 
lar mass distribution of star-forming galaxies, in particular in the Hubble Ultra Deep 
Eield (Eigure[T]). It was then found that the bright patches dominating the optical 
structure were not just random or transient associations of bright stars, but actually 
massive clumps with sizes in the 100-1000 pc range, stellar masses of a few 10^ 
to, in extreme cases, a few 10® M©, and typical stellar ages of a few 10^ yr indicat¬ 
ing relatively young ages (and in particular younger than those of the host galax¬ 
ies) but suggesting lifetimes that are longer than their internal dynamical timescale 
(about 10-20Myr). As random associations would di srupt on such timescales, thes e 
data suggested that these bright regions were bound ( Elmegreen et al.L 20051 2007 ). 
These structures are generally dubbed “giant clumps”, and their host galaxies 
“clumpy galaxies”, although these are just the majority of star-forming galaxies at 
redshift 1-3 and the dumpiness is not a peculiar property of a rare type of galaxy. 
The morphology of clumpy galaxies suggested that these were disk galaxies, based 
in particular on the distribution of axis ratios from face -on to edge-on orientations 
( Elmegreen et ak . 2004; Elmegreen & Elmegreen! 20061) . although the morphology 
of chain galaxies could have been consistent with filamentary alignments of se parate 
small galaxies rather than edge-on clumpy disks ( Taniguchi & Shio\^ 12001 ). 
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Fig. 1 Portion of the Hubble Ultra Deep Field (optical survey) showing more than 10 star-forming 
galaxies at z = 1 — 3 with stellar masses of a few 10*° Mo (largest galaxies on the image). The 
galaxy in the dashed rectangle is a typical example of a Main Sequence galaxy about the mass of the 
Milky Way, located at redshift z— 1-6. Its two apparent neighbours lie at very different redshifts. 
Like many galaxies in this mass and redshift range, it has a very irregular clumpy morphology, 
with a central reddish bulge and a few clumps of 10*“° Mq each. Only the clump to the right of 
the centre is somewhat redder and contains old stellar populations, and only this particular clump 
might be a minor merger of an external galaxy. The others formed in-situ by gravitational instability 
in a gas-rich turbulent disk ^Boumaud et al.Ll2008h and follow a regular rotation pattern around the 
mass centre. The violent instability in the gas-rich medium can also trigger the asymmetry of the 
galaxy through an m = 1 mode, without requiring an external tidal interaction. 
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More robust studies of the nature of clumpy star-forming galaxies were enabled by 
spatially-resolved spectroscopy of the ionised gas, probing the gas kinematics (ve¬ 
locity field, velocity dispersion) as well as che mical abundances an d gradients in the 
interstellar medium. The pioneering study of Genzel et al.l (2006) probed the disk¬ 
like nature of one typical star-forming galaxy at redshift two, with a disk-like veloc¬ 
ity field, and no signature of an on-going or recent merger, in spite of an irregular 
and clumpy morphology. While in a single case a merger might show no observable 
kinematic signature if its velocity field looks like that of a disk, depending on the 
interaction orbit and o n the observer’s line-of-sight, much larger surveys have s ince 
then been assembled ( Forster Schreiber et al. . 20061 20091 : Epinat et al. . 2012 ) and 
quantitative techniques have been used to in terpret the velocity field structure and 


velocity dispersions of the observed systems (iShapiro et al.Ll2008h . These confirmed 


that only a minority of clumpy galaxies display potential signatures from mergers, 
and that clumpier galaxies do not harbour more frequent signatures of recent or 
on-going mergers than smoother galaxies in the same mass range. 

The velocity dispersions ar e high: typically about 50km s~* with la rge spatial 
variations for the Ha gas (e.g., Genzel et al.ll2008t Bournaud et al.lEoO^ . Yet, they 
remain fully consistent with the large scale heights expecte d for these disks if chain 
galax ies, which are relatively thick, are the edge-on version ( Elmegreen & Elmegreenl 


2006h . There are also local irregularities in the velocity field, departing from pure 
rotation with non circular motions up to 20-50km s^*, a few times larger than 
in nearby spirals, but again this is naturally expected from the presence of giant 
clumps: independently of their origin, massive clumps do stir the surrounding ma¬ 
terial by gravitational and hydrodynamic interactions, and they interact with each 
other: models of rotating galactic disks with massive clumps do predict such large 
non-circular motions on kilo parsec scales eve n when the disk is purely rotating 


before clump formation ( Bournaud et al. . 20081) . 


Parametric classifications of galaxy morphology have long been based on low- 
redshift data and models tuned to represent the low-redshift Universe. Yet, such 
classifications now start to be tuned for for high-redshift galaxies, and optimised to 
avoid confusion between internal clumps and mergers. Such morphological classi¬ 
fications confirm the “clumpy disk” nature of the majority of Main Sequence star 
forming galaxies (Eig ure 2) and show close agreement with kinematic classifications 


(ICibinel et al.L 120151 


Signatures of mergers might be more prominent in clumpy ga laxies at intermedi¬ 
ate redshift (z = 0.5 — 1), according for instance to Puechl ( 2010l) . This is a different 
regime in which clumpy galaxies are more rare among star-forming galaxies, most 
of which have already starte d to establish a regular barred spiral structure (in the 


mass range considered here - Sheth et al. 200M Kraliic et al. 2012 ) and have much 


lower gas fractions and densities ( Combes et al.Ll2(^ ). Merger-induced clump for¬ 
mation may then become more prominent compared to higher redshifts. We also 
note that these clumps have lower masses so their dynamical impact may be weaker, 
their response to stellar feedback being likely different, and their role in disk and 
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Fig. 2 Morpholog ical classification of dis ks and mergers at z ~ 2 (Cibinel et al. 201 5). Combining 
the Asymmetry A jConselice et al.Ll2003h and the M 20 parameter jLotz et alil2004h . measured on 
stellar mass maps, is the most efficient way to distinguish clumpy irregular disks and genuine 
mergers. The probability of being a disk or a merger according to these parameters is coded using 
the background colours (note that the scale is logarithmic so disks strongly dominate all colour 
bins, except the last one). The black symbols show a mass-limited sample of star-forming galaxies 
at z 2 in this A — M 20 plane (Cibinel et al. 2015). About two thirds of these objects are secure 
disks, and many others have a high probability of being a disk. The disk fraction is even larger when 
the sample is limited to Ma i n Sequ ence galaxies, excluding the starbursts. Kinematic classifications 
by iForster Schreiber et alj i2009ll are in close agreement with this morphological classification, 
when applied to the same galaxies. Figure courtesy of Anna Cibinel. 


bulge evolution being potentially different as well. For these reasons the relatively 
rare clumpy galaxies below z ~ 1 will not be considered hereafter. 
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2.3 Observational insights on the nature of giant clumps: gas 
content and stellar populations 


The very nature of clumps and their formation process remain uncertain, in spite of 
the fact that their host galaxies are generally rotating disks, with a low frequency of 
mergers. Namely, are these structures formed in-situ, or do they originate from the 
outside in the form of small companion galaxies that have been accreted, or clumps 
of primordial gas that were accreted by the host galaxy before starting to form stars? 

The hypothesis of external star-free gas clouds might be ruled-out by the high 
average density of gas in the cl umps (hundreds of atoms or molecules per cm^, e.g. 


Elmegreen & Elmegreenll^OSh making star formation very efficient, not counting 


the denser s ubstructures likely t o arise given the high observed turbulent velocity 
dispersions dPadoan et all ll999h . 

This hypothesis of clumps coming from the outside as small companion galaxies 
joining a massive galactic disk through dynamical friction requires deep imaging 
to be examined through stellar population studies at the scale of individu al clumps 


( Eorster Schreiber et ail I I 2 OI ll i lWuvts et alll2012 : lElmegreen et al.Ll200^ The vast 
majority of clumps appear younger than expected for small external galaxies at the 
same redshift. Although their young stellar content of giant clumps may bias the 
age estimates by outshining the older stellar populations, the comparison with small 
galaxies at the same redshift shows that the clumps are significantly younger than 
small galaxies. Many clumps have estimated stellar ages of only about 100 Myr with 
no underlying old stellar populations, while such populations would be detectable 
in small galaxies. 

If the clumps really form in-situ in their host galaxy, one should in theory some¬ 
times capture the formation of the clumps during their first internal dynamical 
timescale (<20Myr) and hence some clumps should have extremely young ages 
(about lOMyr only). However, stellar population studies with broad band imaging 


and n o spectroscopy cannot robustly distinguish such very young ages (IWuvts et al 


2012h . Only r ecently such candid ates could be identified with deep imaging and 


spectroscopy ( Zanella et all 2015 ). At the opposite, considering galaxies with stel¬ 
lar masses of a few 10^* ** Mq, a small companion of about 10^ M© should be found 
within a projected distance of 10 kpc for about one third of galaxie|3- If this satellite 
has not been fully disrupted by the galactic tides, its nucleus should be observed 
as a giant clumpQ Such “ex-situ” clumps, with older average stellar ages and an 
underlying old population, are indeed found in some cases. E or instance the rep¬ 
resentative clumpy galaxy dissected in iBournaud et al.l (120081) contains one clump 
which is much redder and older than the others, and also exhibits a larger devia¬ 
tion from the underlying disk velocity field, making external origin most likely for 


* Tills estimate is simply based on the mass function of galaxies and assuming a random geomet¬ 
rical distribution of satellites within the virial radius. 

^ With a kinematics that could become preferentially consistent with that of the host galaxy disk 
through gravity torques and dynamical friction within one galactic dynamical time, i.e. about 
100 Myr. 
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this one. Other candidates are foun d in iForster Schreiber et al.l (120111) and appear in 
the statistics of IWuvts et al.l (l2012r) . Nevertheless, such ex-situ clumps remain rela- 
tively rare: most observed clumps are actually different, with youn ger stellar ages 
( Elmegreen et ah . 20091: Bournaud et all 2008 : Wuyts et all 2012 ). The ability to 
identify such ex-situ clumps is actually reassuring that the non-detection of such old 
populations in the other clumps is robust, and probes their recent in-situ formation. 

Since most of the clumps formed recently inside their host galaxy, and giyen that 
they are gravitationally bound (based on the observed stellar masses and velocity 
dispersions, e.g., Elmegreen et al. 2005, Genzel et al. 2008), their formation likely 
involves a gravitational (Jeans) instability in a rotating disk, sometimes also called 
Toomre instability. Eor such instabilities leading t o the formation of bound objects 
in a rotating disk to arise, the key requirement is a Toomrel ( 1964t) parameter Q be¬ 
low unit}0. Based on the observed rotation velocities and velocity dispersions (see 
previous sections) this typically requires gas density of the order of IOOMqPC^^, 
an order of magnitude larger than in nearby disk galaxies, and implying interstellar 
gas masses comparable to the stellar masses, i.e. gas mass fractions of about 50% of 
the baryonic mass. Such high gas densities had long been foun d only in the stronges t 


starbursts galaxies (likely merger-induced) at any redshift (e.g.- lTacconi et alll2008l) 


but not in normally star-forming galaxies. The discovery that Main Sequence galax¬ 
ies at redshift z > 1 are actually v ery gas-rich with gas fractions of abou t 50%, 
just c ounting the molecular content ( Daddi et al. . 2008 . 2010l: Tacconi et al. . 20ld 
20131), has shed a new light on this issue. These high gas fractions are estimated from 
CO line observations, and hence subject to uncertainties on the conversion of CO 
luminosity t o H? mass. Yet, new data probing the CO molecule spectral line energy 
distribution ( Daddi et all 2015 ). compared t o detailed modelling o f the CO exci¬ 


tation and emission in high-redshift galaxies (IBournaud et al.L 12015h . confirm high 


luminosity-to-mass conversion ratios and high gas mass fractions. Eurthermore, the 
high g as fractions are also confirmed by independent estimates based on dust prop¬ 
erties ( Sargent et al. . 2014; Magnelli et al.L 2014; Genzel et al. . 2014). 

Hence, the high inferred gas surface densities lead to Toomre parameters around 
unity for velocity dispersions of 30-50 km s^\ and the Toomre Q parameter would 
be even lower if the velocity dispersions were only a few kms^^ as is the case in 
nearby spirals. Thus, axisymmetric gravitational instability should arise. All numer¬ 
ical experiments modelling disks with masses, sizes and gas fractions representa¬ 
tive for the high-redshift star-forming galaxi es consider here do show clump for¬ 
mation through gravit ational instability (see Noguchil 1999t Immeli et al. . 2004a : 
Bournaud et al.L 2007 . and the more detailed models discussed hereafter). The grav¬ 
itational stirring of the gas ensures that the turbulent velocity dispersions do not 
stay below the observed level of ~50kms^*: namely, the gas turbulent motions 
may also be powered by infall and stellar feedback, as we will review in the next 
sections, but at least the release of gravitational energy through the instability is suf- 


^ Although the Toomre Q parameter is strictly meaningful only in an axisymmetric disk before 
strong perturbations arise. No te also that in a thick disk the instability limit is about 0.7 rather than 
Q < I IBehrendt et ali[2014l and references therein). 
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ficient to maintain high turbulent dispersions and self-regulate the disk at a Toomre 
parameter 2 ~ 1. 

The high velocity dispersions imply that the Jeans mass (or Toomre mass) is 
high, typically 10^ — 10® M©. This sets the high mass of the giant clumps form¬ 
ing through the associated instability. The prop erties of the i nstabili ty in these disks 
was also modelled through analytic models by Dekel et ^ ( 2009bll . with results in 
close agreement to that of numerical models, and simpler analytic estimates. The 
giant clumps observed in high-redshift galaxies can thus be considered as the di¬ 
rect outcome and signature of the disk instability with a high characteristic mass, 
robustly expected from the basic observed properties (mass, size, rotation speed 
and gas fraction) of these high-redshift galaxies. This probably applies to most of 
the observed giant clumps, except the few ex-situ clump (minor merger) candidates 
with older stellar ages (see above). Note that the gravitational instability actually 
arises in a two-component disk with r oug hly half of its mas s in gas, and half in 
stars: we refer the reader to Jod ( 19961) and Elmegreen ( 201 ih for theoretical work 
on the two-component stability, and Behrendt et alj( |20l4' for a detailed analysis of 
the disk stability in numerical simulations. The process is qualitatively unchanged 
compared to the 2 ~ 1 self-regulated instability in a single-component disk de¬ 
scribed above; furthermore the gaseous and stellar velocity dispersions are probably 
about similar in these young galaxies (e.g., Bournaud et al. 2007) in which case the 
single-component stability analysis directly applies. 


2.4 The formation of gas-rich clumpy unstable galaxies in the 
cosmological context 


The irregular and clumpy structure of high-redshift star-forming galaxies is the out¬ 
come of their high gas fractions and densities. Theoretically, these high gas fractions 
are explained by the high rates of external mass infall, which are not compensated 
by excessively high star formation rates consuming the gas reservoirs, but rather 
preserved in a long-lasting steady state for Main Sequence galaxies. Recent cosmo¬ 
logical models have highlighted the fact that high-redshift galaxies mostly accrete 
their baryons in the form of cold diffuse gas rather tha n hot reservoirs and compan¬ 
ions galaxies ( Dekel et al. . 2009al : Brooks et al. . 2009h . which further helps the gas 
to rapidly join the cold star-forming disk. At the opposite, a too large contribution 
of galaxy mergers in the cosmological galaxy growth budget would form a massive 
stellar spheroid (bulge or stellar halo) too early and would stabilize z — 2 galaxies 
against giant clump formation, with a much low er turbulent speed and character istic 
mass for any residual disk instability at Q ~ 1 ( Bournaud & Elmegreeni 2004 . In 
detail, the cold accretion streams do not directly join the star fo rming disk (see Fig- 
ure they might be affected by the hot circum-galactic gas ((Nelson et all 1201 3h 
but more importantly need to dissipate the high kinetic energy from the infall, possi¬ 
bly by turbulent dissipatio n in circum-galactic regions ( Elmegreen & Burkerl , l2010t 
Gabor & Bournaud 12014 ) and/or forming extended rotating reseryoirs that gradu- 
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ally feed the dense star-forming disk (iDanovich et al.L 120141) . Yet, detailed cosmo 


logical simulations show that part of the streams can direc tly feed the central few 
kiloparsecs of galactic disks in less than one dynamical time (IDanovich et al.L 120141) . 

As a result of the high accretion rates, galaxies in cosmological simulations 
around redshift two have high gas fractions, although modern simulations still have 
difficulties to preser ve sufficient gas reservoirs b y avoiding excessive star formation 
at early epochs (see Dekel & Mandelkei . 20l4 and references therein). Neverthe¬ 
less simulations with detailed stellar feedback models do produce some steadily 
star-forming galaxies with up to ~ 30-40% of gas at z = 2 — These simulations, 
further detailed in the next Section, display the expected gravitational instability for 
such gas-rich disks, and produce giant clumps and other dense features by gravi - 
tational instability in these disks ( Agertz et ah . 20091: Ceverino et al.L 201o[ 2012). 
Only a limited fraction of the clumps are “ex-situ” clumps, resulting from the accre - 
tion of small nucleated companions or external gas clumps (IMandelker et al.ll2014h . 
fully consistent with the observations reviewed above. 


Fig. 3 Primordial galaxy 
fed by three cold streams 
of gas, in an idealised high- 
resolution simulation based 
on typical parameters mea¬ 
sured i n cosmological simu ¬ 
lation s l lGabor & Boumaud . 
|2Q14) . The flows join the disk 
through a turbulent interface 
with extended circum-galactic 
reservoirs, before the gas dis¬ 
sipates its energy and feeds 
the cold star-forming disk, 
which keeps a high gas frac¬ 
tion, high velocity dispersions 
self-regulated at a Toomre pa¬ 
rameter (2 — 1, and a clumpy 
irregular morphology (image 
size: 100x70 kpc). 



3 Mechanisms of bulge growth through high-redshift disk 
instabilities 

Knowing the properties of high-redshift star-forming galaxies from the previous 
Section, in particular the fact that they are subject to a violent clump instability 
rather than a weak axisymmetric instability as in nearby barred spirals, we now 

^ Although these total gas fractions seem to remain lower than the observed molecular gas frac¬ 
tions, especially if these are the most gas-rich galaxies in simulated samples. 
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review the mechanisms through which bulge growth can be triggered and regulated 
by this instability mode. 


3.1 Clump migration and coalescence 


With a mass of a few 10^^® Mq, a giant clump in a high-redshift disk galaxy could 
behave just like a dwarf companion galaxy of similar mass, except that being dark 
matter free makes its mass distribution more spatially concentrated. It will in partic¬ 
ular undergo dynamical friction on the underlying gaseous and stellar disk, and dark 
matter halo, and through this process will dissipate energy and angular momentum 
through increasing the velocity dispersion (i.e. the internal kinetic energy) of these 
massive components. In a rotating disk this will lead to inward migration of each 
giant clump until it reaches the galaxy centre - like in a minor galaxy merger. In ad¬ 
dition, as they lie in the disk plane, clumps also undergo gravity torques from other 
regions of the disk. A clump that forms in a purely rotating disk will break the sym¬ 
metry of the mass distribution and induce a kinematic response in the form of a spira l 
arm or tidal arm denser than the average disk, in the disk plane ( Bournaud . 201(lh . 
This over-dense region will then exchange angular momentum with the clump it¬ 
self. If most of the mass lies at radii larger than the clump in the galactic disk, as 
is the case as soon as the disk is sufficiently extended radially, the strongest gravity 
torques will be from this arm onto the clump. Given that the outer disk has a slower 
angular velocity than the clump, this arm is trailing with respect to the rotation of 
the disk, so that the gravity torques exerted on the clump are negative, accelerating 
the inward migration already resulting from dynamical friction. This torquing pro¬ 
cess is most efficient in a gas-rich disk as the cold gas component makes the tidal 
arm response stronger than in a pure stellar disk. 

As the process involves both dynamical friction and gravity torques, many nu¬ 
merical simulations have been used to stu dy clump migration and estim ate the mi- 
gr ation timescale, starting with those of IShlosman & Noguchil (Il993h and those 
of Noguchi ( I 999 I) . the latter being directly motivated by the first observations of 
chain galaxies and suggestions that the latter were internally fragmented disks by 
Cowie et al. ( 19961) . 

A more de tailed treatment of the hydrodynamics and interstellar gas physics was 
introduced by Immeli et al. ( 2004alf9) . and the simulations of Bournaud, Elmegreen 
&. Elmegreen (2007) were designed to finely correspond to the observed proper¬ 
ties of star-forming galaxies in the Hubble Ultra Deep field at redshift z = 1 —2 
with stellar masses of lO*** to IO'^Mq (see Eigure|4]i. In the latter, the migration 
timescale of giant clumps to the galaxy centre is of the order of 300-500Myr, with a 
large scatter depending on the clump initial formation radius, and also of its interac¬ 
tion with other giant clumps and dense features in the disk. All of the experiments 
cited above are idealised models of isolated galaxies, lacking external replenish¬ 
ment of the disk by cosmological gas infall, accretion of smaller galaxies, and a few 
bigger mergers. As a result of star formation the gas fraction gradually decreases. 
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Fig. 4 Simulations of a gas-rich (50% gas fraction) disk galaxy with initial parameters represen¬ 
tative for star-forming galaxies at z ~ 2. The unstable gaseous disk forms a ring that quickly 
fragment into giant clumps. The clumps migrate inward and coalesce into a bulge, while the stel¬ 
lar disk is significantly thickened, and the disk radial profile, initially flat, is re-distributed into an 
exponential. The bulge formed here is a classical bulge with a Sersic index of 3.5^.0. Blue codes 
gas-dominated r egions and red codes s tar-dominated ones. Snapshots are separated by 100 Myr. 
Simulation from iBoumaud et^ i2007h . 


and this c an lead to over-estima ting the clump migration timescale, as noted for in¬ 
stance by Ceverino et al. ( 201(]ll . However the clump migration process is so rapid 
(for clump masses above 10 ^ Mq at least) that the gas fraction decreases by less 
than a third of its initial value over this period, which is not larger than other un¬ 
certainties like the fact of not being able to observed any reservoir of atomic gas 
in these galaxies. Indeed cosmological simulations with external mass infall have 
reproduced the clump formation and migration process and find timescales that are 
consistent with the above, or only slightly shorte:[§. 


^ Note that shorter migration timescales in cosmological simulations may also arise if the galaxies 
are too compact, or have too concentrated dark matter halos enhancing the dynamical friction 
process. Actually, the cosmological simulations of redshift two galaxies tend to have too low gas 
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An observational signature of inward clump migration, if they survive stellar feed¬ 
back (see below), could be an age gradient with, on average, older cl umps found at 


smalle r radii. This is extensively quantified in the simulation sample of lMandelker et al 


(l2014ll . In detail clump migration does not imply that young clumps cannot be 
found at small radii: external gas can feed high gas fraction in the inner disk 
( Danovich et all 20141) and new clumps can form at small radii in gas-rich disks, 
except in the central 1-2 kpc that are stabilis ed by strong shear and bulge mass 
( Bournaud et all 2007 : Mandelker et all 2014). The expected signature is rather an 
absence of aged clumps in the outermost disk, and clumps formed there have mi¬ 
grated inward in 100-200 Myr (exceptions could still be found for moderate-mass 
clumps that can be scattered out to large radii by bigger clumps). 

Observationally, statistical samples or resolved clumps remain limited, and their 
ages are hard to estimate. Not only the age of stars in a clump is not a direct tracer 
of its age (because clumps loose and re-accrete material, see next Section), but also, 
the stellar age estimators are strongly dependent on many pa rameters such as the as - 
sumed star formation histories, especially at high redshift ( MMUston_etaLL201^^ 
Never theless, an age gradient is tentatively observed by iForster Schreiber et al 


(1201 Ih . in quantitative agreement with clump migration and central coalescence 
within a timescale of at most 500Myr (see also lGuo et al.Ll2()lZl2014h . 


3.3 Stellar feedback, outflows, and the clump survival issue 


A key issue in the process of inward clump migration (and subsequent coales¬ 
cence into a central bulge) is their response to stellar feedback. In nearby galax¬ 
ies molecular clouds are estimated to have short lifetimes of the order of 10- 
20 Myr, under the effects of supernovae explos ions and other feedback processes 
( Hennebelle & Falgarone . 2012 : Murravi 2011 ). While giant clumps are typically 
a thousand times more massive than the Milky Way biggest gas clouds, they are 
also ten times larger in each dimension, so their 3-D mass density is not neces- 
sarily much high e r. As t hey form stars at high r ates of a few Mf^ vr~* pe r clump 
( Elmegreen et al. . 2007t Forster Schreiber et al. . 2009t Wuvts et aH 2012) the re¬ 
leased energy per unit gas mass is of the same order as that found in nearby star¬ 
forming clouds, or slightly higher. This raises the important question of clump sur¬ 
vival against stellar feedback. In particular, having a released feedback energy per 
unit gas mass of the same order as in nearby molecular clouds does not mean that 
the giant clumps will be disrupted in a similar way or on a similar timescale: their 
gas also lies in a deeper gravitational potential well. 


fraction because of some largely unexplained early consumption of the gas (e.g., Ceverino et al. 
2012, Keres et al. 2012, with typical gas fractions at best around 30% at z = 2 even counting all 
the cold gas within a large radius). 
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A first attempt to address this issue in Elmegreen et alJ (2008) concluded that 
if feedback was strong enough to disrupt the giant clumps within their migration 
timescale, it would also severely thicken the gas disk and heat the stellar disk well 
above the observed levels, almost disrupting any pre-existing rotation-dominated 
stellar disk. This was however based only on energetic supernovae feedback, while 
other stellar feedback mechanisms may be more likely to disrupt clumps in a realis¬ 
tic way, without disrupting the whole host galaxies. In particular, radiation pressure 
from young massive stars on the surrounding gas and dust may inj ect enough mo¬ 
mentum into the clumps to disrupt the clumps ( Murray et all 201 ll) . 

It has long remained difficult to address this issue in numerical simulations, 
mostly because stellar feedback can only be modelled through uncertain sub-grid 
models, even if modern hydrodynamic simulations of galaxies can reach sub- 
parsec spatial resolu tions with mass resolution elements of the order of IOOMq 


(iRenaud et al.L 120131) . In fact, even the star formation rate which determines the 


powering rate of feedback relies on sub-grid models. Even if the star formation rate 
of entire galaxies or giant clumps is realistic compared to observations, changing the 
sub-grid model may significantly alter the spatial distribution of star formation, es¬ 
pecially in resolution-limited simulations. A reassuring point is that state-of-the-art 
idealised simulations of galactic physics can now robustly model gaseous structures 
up to densities of 10® cm^^ or more without being at their spatial resolution limit yet 
(i.e., the Jeans length at these densities can remain larger than a couple of resolution 
elements without necessarily requiring to bias the hydrodynamics by adding a tem¬ 
perature or pressure floor, Renaud et a l. 2013). The fact that stars form with a quasi- 
universal efficiency in such dense gas ( Krumholz & McK^ 2005; Gao & Solomon . 


2004tlGarcfa-Burillo et al.Ll2012l) implies that at least the first step of star formation. 


which is dense gas formation, is explicitly resolved in these simulations; the subse¬ 
quent sub-grid modelling of star formation at fixed efficiency in high-density gas is 
consistent with observations down to scales much smaller than that of giant clumps 
- this is now achieved in idealised simulations, but unfortunately remains out of 
reach of cosmological simulations so far. 

The modelling of supernovae feedback is highly uncertain, in particular because 
it is often modelled through thermal dumps of the released energy heating the sur¬ 
rounding gas, while real supernovae remnants include a large frac tion of the en¬ 


ergy i n non-thermal processes, which dissipate on slower timescales (ITevssier et al 


2013h . Eurthermore, models including the other kind of feedback processes such 


as 


stellar winds, photo-io nization, and most importantly radiation pressure were de¬ 
veloped only recently (IHopkins et al.l 120131: iRenaud et al.l 1201 3l see Eigure |5]l and 
include free parameters, an important one for radiation pressure being the number 
of scatteri ng events that a pho ton can undergo in gas cloud before escaping from 


the cloud (iMurrav et al.l 1201 Ih . Another key parameter is the initial mass loading. 


namely whether the available energy or momentum is diluted into a large or a small 
mass (and volume) of gas. This parameter has long been unresolved in numerical 
simulations, an d sometimes adjusted to generate ad hoc galactic outflows and study 
their fate (e.g., Oppenheimer & Dave . 2006t Genel et ^ 2012h . The highest reso¬ 
lution simulations of galaxies now become capable of resolving the typical distance 
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Supernova Hll+radiation pressure all feedback 


Fig. 5 Edge-on views of three simulations of the same gas-rich clumpy galaxy that has been 
evolved with different stellar feedback models during the last 80 Myr (from left ro right: supemovae 
only, photo-ionization and radiation pressure only, and all mechanisms together, respectively). The 
gas density is shown and the outflow rates are indicated in the panels in Mq yr“ * (measured 2 kpc 
above/below the disk mid-plane). Outflows are launched by the giant clumps, and the models show 
the strongly non-linear coupling of feedback mechanisms: the total outflow rate in the simulations 
with all feedback processes together is well above the sum of the outflow rate in the independent 
cases. Similar non-linear coupling was noted by Hopkins et al. (2013). These simulations use the 
feedback models proposed by Renaud et al. (2013) and are similar to those presented in Boumaud 
et al. (2014), with 3 pc spatial resolution. 


over which photons from young stars redistribute their momentum into the ISM 
and start to estimate this loading factor from physical principles, even though ex¬ 
plicit radiative transfer calculations robustly resolving these typical scale lengths are 
out of reach from galaxy-scale models and become feasible only in cloud-scale or 
clump-scale simulations. 

Some models of gas-rich galaxies with intense feedback have found that the giant 
clumps could be short lived, even with clump masses of the or der of 10^ M^. Thi s 
is the case for instance in the cosmological sim ulations from iGenel et al.l (l2012h . 
or the idealised models of Hopkins et al.l ( 2012 ). It is nevertheless remarkable that 
in these short-lived clumps models the clump lifetimes are very short, not larger 
than 50 Myr, hence appearing inconsistent with the stellar ages estimated for real 
clumps, often reaching 100-200Myr and more (see above and Wuyts et al. 2012). 
In these models, the clump disruption is obtained in one or two generation of star 
formation and evolution, rather than through gradual, steady outflows on the longer 
term. Models with strong feedback and no long-lived clumps actually tend to lack 
giant clumps at all, strongly reducing and mass and/or number of clumps formed, 
at the point of being inconsistent with forming the majority of observed clumps 
by in- situ instability, as highlighted recently in the simulations of Tamburello et al.l 
(I 2 OI 4 I) . However, in such models where in-situ clump formation is suppressed, a 
different (ex-situ) origin of clumps is not explained: in particular their stellar pop¬ 
ulation ages can hardly be reconciled with minor mergers - minor mergers can ac¬ 
tually b e identihed as a sub-population of clumps tha t contain older stellar pop¬ 
ulation ( Bournaud et al.l 2008 : Elmegreen et al. . 2009) and these are only a small 
fraction of giant clumps. The suppression of in-situ giant clump formation obtained 
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in the models of iTamburello et alJ (120141) could in fact result of the low surface den¬ 
sity of the disks in their initial conditions, which were inspired from cosmological 
simulations (which in turn may consume the disk gas too early), rather than from 
observ ed gas surface densi t ies estimated throug h detailed analysis of the dust prop¬ 


erties (ISargent et all 12014t iGenzel et all 120141) and carbon monoxide spectral line 


energy distribution studies ( Daddi et al. . 2015h Hence, a common drawback of all 
theoretical models without long-lived clumps, is that either clump formation is sup¬ 
pressed, or the clump formation/disruption cycle is very short (<50Myr): in any 
case this appears inconsistent with the observations that commonly probe clump 
stellar ages of 100-200Myr and more. On the other hand, for long-lived clumps in 
a typical star-forming galaxy of stellar mass lO**'^" M©, the migration timescale 
from the clump birth site to the galaxy central kpc should be 300-500 Myr, which 
appears slightly longer than the observed average stellar ages in giant clumps: this 
led Wuyts et al. (2012) to mention that clump disruption might be faster than clump 
inward migration. Yet, the clump stellar population ages provide only a lower limit 
to the clump real ages (see next paragraphs in this Section). 

Actually, simulations with a thorough accounting of stellar feedback processes, 
including not just supernovae but also radiation pressure and other feed back mecha 


nisms, do n ot necessarily predict short-lived clumps. In c o ntrast with iGenel et al 


(201J) and Hopkins et al. ( 201^. models in Ferret et al. ( 2014 ). Bournaud et al 


(120141) or ICeverino et al.l(l2014j) include non-thermal and radiative feedback schemes 
and do find long-lived giant clumps - at the same time they do correctly predict short 
lifetimes for gas clouds below 10^ Mf^i like in low-red shift galaxies. A different ap¬ 
proach to feedback modelling by Perez et al. ( 2013h also find long-lived clumps 
for any acceptable amount of stellar feedback, and even when strong outflows are 
launched by the giant clumps and their host galaxies, w ith o utflow rates consiste nt 
with observations such as those of Newman et al.l ( 20121) and Genzel et al. ( 2011 ). 


Important constraints on the lifetime of giant clumps lifetimes and their abil¬ 
ity to migrate inward toward bulges result from the fact that giant clumps are not 
quasi-closed-box entities, but steadily exchange mass with the surrounding medium 
in the host galaxy, through outflows and inflows of both gas and stars. Hence the 
ages of stars inside a given clump are not equal to the age of that clump since its 
actual formation: clumps have a sort of wave-like behaviour, although the pattern 
speed of the m = 0 instability is almost equal to the disk rotation speed. Clumps 
may loose gas through stellar feedback, but more generally loose material through 
gravitational tides. At the clump half-mass radius, the gravitational force from the 
entire clump is only a few times larger than that from the entire galaxy. In other 
words, clump densities are only marginally higher than the limiting tidal density 
( Elmegreen & Elmegreenl 12005 ) and clumps grad ually loose aged stars b y dynami¬ 
cal eyaporation toward the galactic potential well ( Bournaud et al.L 12007 ). 

Conyersely, the clumps haye a large cross section (of the order of O.l-l.Okpc^), 
and wander in a disk that contains substantial amounts of gas, eyen outside the giant 
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clumps themselvefl Given this large cross-section of clumps, their relative veloc¬ 
ity of 10-50kms^* with respect to surrounding gas, and a density of ^10cm^^ 
for the latter, accretion rates of I-IOMq yr^* onto each giant clump are expected 
for pure ballistic capture. The gravitational potential well associated to the giant 
clumps may actually e nhance the accretion . The f irst detailed estimates of this pro¬ 
cess were provided by Dekel&jCrumholz (2o3. Detailed hydrodynamic simula¬ 
tions with an accurate AMR code (lTevssierll2002ll . very high resolutions of 3-6 pc, 
and detailed feedback models c ombining supernovae, photo-io nization and radia - 
tion pressure, were presented in Bournaud et al.l ( 2014 . see also Ferret et ^12014 ). 
These simulations confirmed that, independently on the details of stellar feedback 
and its “strength”, clumps accrete a few solar masses per year of fresh gas. This 
gas accretion onto the clumps roughly compensates for both the gas consumption 
through star formation and the losses by gaseous outflows and dynamical evapo¬ 
ration of aged stars. This means that the clump actually evolves in a steady state, 
which can be d escribed by a so-call ed “bathtub model” more commonly used for 
entire galaxies ( Bouche et al. . 2010 h : the gas infall rate is equal to the sum of the 
star formation rate and gas outflow rate, keeping the total mass constant, letting the 
system evolve in a steady state. The idea behind the steady state regulation is that 
any increase in the gas infall rate will be compensated for by the star formation rate, 
which has a non-linear response, and vice-versa for any decrease in the gas infall 
rate. The clump mass is then about stabilised, with some fluctuations, around its 
“initial mass”, i.e. its baryonic mass after 1-2 times its own gravitational free-fall 
time (10-20Myr). 

An important prediction of the long-lived clump scenarios is that the average 
age of stars contained by a given giant clumps is younger than the actual clump 
age, measured since its formation by gravitational collapse of the gas-rich disk. The 
clumps experience a moderate starburst during thei r first 10-20 Myr bef ore feedback 


regulates star formation in a steady state regime (IZanella et al.L 120151) . Then, they 


continue to form stars steadily at a higher rate than a closed-box system, due to 
(re-)accretion of gas from the larger-scale galactic reservoirs. This keeps the average 
stellar age younger than the clump age. Furthermore, aged stars leave the clump 
gradually throughout the effects of dynamical heating and evaporation toward the 
galactic tidal held. As a result the stellar age is even younger than the clump age, 
and typical stellar ages of 100-200 Myr are predicted for clumps of real age 300- 
500 Myr, as the age of the stellar constant of giant clumps tends to saturate at 200- 
250 Myr (Bournaud et al., 2014). 


® The presence of large amounts of gas between the giant clumps cannot be mapped spatially in 
CO surveys yet, but is predicted in the idealised and cosmological simulations of gas-rich un¬ 
stable disks cited above, and confirmed by two observational arguments: (1) the emission from 
young stars in the u ltraviolet contains a wide s pread component behind the giant clumps, tracing 
relatively dense gas lElmegreen & Elmegreenl [20051 ) and (2) the CO spectral line energy distribu¬ 
tion has two components, a high-excitation one attributable to dense clumps, and a low-excitation 
one corresponding to lower-density, large-scale background gas reservoirs iDaddi et all l20ll 
iBoumaud et alil2015l) . 































18 


F. Boumaud 


If clumps were disrupted by stellar feedback-driven outflows, the same processes 
of mass loss and accretion onto the clumps, which are driven by gravitational dy¬ 
namics, would still be present. Hence the stellar ages would still set a firm minor 
limit to clump ages. Hence the observed stellar ages of g iant clumps, typically of 
at least 100-200Myr ( Wuvts et al. . 2012 : Guo et all 2012 ). constitute evidence that 
clumps are not disrupted in a few tens of Myr as predicted by some feedback models 
- actually, all the short-lived clump models reviewed above predict short lifetimes 
smaller than 50 Myr. The observed stellar ages appear only consistent with models 
of clump survival and migration toward the galactic center, the ultimate limiting 
factor to the clump lifetime being their coalescence with other clumps or with the 
galactic bulge, explaining that Gyr-old giant clumps are not observed either. 


A typical star-forming galaxy at redshift z = 1 — 3 can thus be expected to experi¬ 
ence the migration and central coalescence of giant clumps of 10 *^® Mq of gas and 
starfl Given the observed (and simulated) number of clump per galaxy, combined 
with their theoretical lifetimes and observed stellar ages, the central coalescence of 
a giant clump should typically occur at a rate of lOGyr^* for a galaxy of stellar 
mass 10**'^** Mq, i.e. one giant clump every 10^ yr or so. If the unstable steady 
state lasts 2 Gyr, this means that the baryonic mass reaching the bulge can be of 
order of lO*** Mq - or even higher without strong stellar feedback, as in such a case 
the giant clumps have masses that increase through the accretion of surrounding gas 
without outflow regulation. 


While the mass reaching the bulge can seem very high, we will see later that it 
does not necessarily entirely fuel a too massive bulge: in particular, a large fraction 
of the mass is still gaseous and this component can be expeled outward and/or form 
a central rotating disk rather than a bulge. In the next parts we first examine below 
the structural properties of the bulges formed by the central coalescence of giant 
clumps, and will review the bulge mass fraction issue subsequently. 


3.4 Instability-driven inflows 

Another mechanism associated to giant clumps and disk instability, but different 
from giant clump migration and coalescence itself, can also grow the central mass 
concentration of high-redshift disk galaxies, and potentially the bulge mass. Gi¬ 
ant clump and other dense features formed by gravitational instability exert gravity 
torques on the rest of the disk’s gas, which transfers angular momentum. Clumps 
are located close to their own corotation radius (or slightly inside their corotation if 
dynamical friction has slowed down their rotation speed compared to the rest of the 
disk material). Material located at smaller radii in the disk thus rotates faster than 
the clumps, in terms of angular velocity. It then responds mostly as a leading tidal 
arm, found on the leading side of the clump compared to the galactic rotation. The 

^ Note that the clump.s remain gas-rich as they re-accrete gas and lose aged stars, which compen¬ 
sates for the gas depletion through star formation and gaseous outflows. 
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presence of multiple clumps and other features can make the pattern hard to identify 
but a striking example with respect to one of the main clumps in a clumpy galaxy 
simulation is shown in Figure |6] 



Fig. 6 Face-on view of the gas in a high-redshift galaxy simulation (image size: 8x 12kpc). The 
galactic rotation is counter-clockwise. Note the spiral armlets which are often on the leading side 
of giant clumps inside the clump radius, and on the trailing side in the outer disk. Gravitational 
torques from the clumps onto this inter-clump gas drive a continuous inflow of gas toward the 
galaxy center. 
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The material on the leading side of giant clumps undergoes negative gravity 
torques and loses angular momentum. The material in the outer disk gains angu¬ 
lar momentum in ex change. The process is similar for any instability that breaks the 
disk symmetry (e.g., Combes & Gerin . 1985 : Bournaud et al.L 2005 ) but the gravity 
torques in the case of clump instabilities at high redshift are typically 10-20 times 
larger than for low-redshift secular instabilities (namely, spiral arms and bars): it can 
even transfer outward 100% of the initial angular momentum in just one rotation pe¬ 
riod ( Bournaud et al. . 201 ih. compared to 5-1 0% per rotation period for strong bars 
in low redshift spirals (iBournaud et al. . 2005h . The corresponding mass inflow rate 


or 


for a typical high-redshift star-forming galaxy is then of the order of 10 M© yr 
more. 

The gravity torques between clumps and non-axisymmetric features are the main 
mechanism through which gravitational energy is pumped into the interstellar gas. 
As reviewed in the previous Section, these unstable disks evolve in a self-regulated 
regime at 2 ~ 1 with high velocity dispersions <7 50km s^*. Turbulent en¬ 
ergy in the interstellar medium typically dissipates in a local crossing-time and 
more energy needs to be pumped in the turbu lent cascade for steady-state systems 
(e.g.. lMac Low . 1999 : Bournaud et al.L 2011). The specific energy loss rate is then 
(7^/ (2t) where T is about 10 Myr, the energy being dissipated mainly through small- 
scale compression and shocks that heat the gas, which subsequently radiates the 
energy away. The radiative losses are balanced by the global inflow of gas down 
the galactic gravitational potential at an mass inflow rate M, releasing an energy 
rate MV^/2 for a galactic circular velocity Vc- Typical outflow rates estimated for 
high-redshift star-forming galaxies are of the order of ~ IO M^ti vr~^ to compensate 
for the turbulent dissipation in a steady state: = MV^ (lElmegreen & Burkert . 


tor me turbulent aissipation in a steaay state: o 
2 OI 0 I: Bournaud et all 2011 : Genel et al. . 2012 ). 


Studies of the instability-driven inflow (see for insta nce iKrumholz & Burkerl , 
2 OI 0 I: Elmegreen & Burkert , 2010 l: Bournaud et al. , 2011) highlight the fact that the 
inflowing gas is not fully consumed by star formation. Actually, a large fraction of 
the initial inflow, several solar masses per year, typically flow onto the central kpc 
region or “bulge region”. The gaseous inflow will not necessarily feed a classical 
bulge as the material is dissipative and discy, but feeds central star formation at a rate 
of a few solar masses per year. This is consistent with observations showing that this 
“bulge region” (about the central kpc) has younger stellar populations and more sus¬ 
tained star formation in the most unstable/clumpy galaxies than in smoother disks 
of similar mass and redshift ( Elmegreen et al.Ll2009[ l2(^ ). The instability-driven 
inflow thus increases the central concentration of gas and young stars, this stellar 
mass being potentially scattered into a pressure-supported spheroidal bulge during 
subsequent relaxation events (which can include: coalescence of other giant clumps, 
major interactions or minor mergers). Eew simulations have studied in the outcome 
of the central mass concentration grown by this inflow, compared to direct bulge 
growth by central coalescence of clumps, mostly because the two processes would 
be hard to distinguish. Yet, it appears clearly that the stellar mass gathered in the 
central kiloparsec by the clump migration and instability-driven inflows ends-up in a 
bulge-like structure rather than just the innermost regions of a radially-concentrated 
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rotating disk: it comes in excess of the disk exponential m ass profile, and has high 


velocity dispersions dominating over the res idual rotation (lElmegreen et all 12008 


Inoue & Saitoh , 2011 : Bournaud et ah . 201 ih . 


The absence of violent relaxation in the process of instability-driven gaseous 
inflows should produce only a so-called pseudo-bulge, namely a low Sersic in¬ 
dex structure with substantial residual rotation. Yet, subsequent relaxation through 
clumps and mergers can make this mass contribute to a so-called classical bulge, 
i.e. a highly concentrated structure with virtually no angular momentum left. This 
global picture has not been studied in detail for disk-dominated galaxies with the 
mass of the Milky Way or up to 10^' M©, but it has been studied for more mas¬ 
sive galaxies. As detailed in Section 5, these massive galaxies become compact 
spheroid-dominated throughout the violent disk instability process. The global in¬ 
stability driven inflow plays a major role in turning the initial disks into compact 
concentrated objects, while relaxation induc ed by giant clumps and some mergers 
turn them into “classical” spheroids (see also Zolotov et al.ll2014 ). 


3.5 Properties of bulges from high-redshift disk instability 


The properties of bulges resulting from the violent instability of high-redshift disk 
galaxies remain uncertain, as they largely depend on the lifetime and evolution of 
clumps against stellar feedback. If the clumps are short-lived, disrupted by feedback 
faster than their inward migration timescale, there is still a diffuse inflow of inter¬ 
clump gas driven by the instability (see above, Hopkins et al. 2012, and Bournaud 
et al. 2011), which will grow a low-concentration pseudo-bulge if not another relax¬ 
ation process affects the central region. At the opposite, the models with long-lived 
clumps, more consistent with the observed ages and outflow rates. And that the insta¬ 
bility can grow classical bulges, wi th a low rotationa l suppo rt and high Sersic indices 
after repeat ed clump coa l escence ( Elmegreenet al.L 2008h and short star-formation 
timescales ( Immeli et al.L 2004bh . Bournaud et al. (12007 ) have shown that for mod¬ 
els scaled to the observed properties of galaxies in the Hubble Ultra Deep Eield, the 
disk material is redistributed into an exponential profile during clump migration and 
classical bulge growth. 

Nevertheless, these first models included only supernova-like feedback schemes 
and lacked a more complete accounting of stellar feedback processes. The full se¬ 
ries of stellar feedback process es regulates the mass and gas richness of clumps 
even if they remain long-lived. Inoue & Saitohl ( 2012h proposed that the resulting 
bulge could rather be a pseudo-bulge, even if it is an old and metal-rich one. It 
nevertheless still seems possible to grow highly-concentrated classical spheroids, 
with detailed and thorough stellar feedback models, at least at high galactic masses 
( Ceverino et al. . 20141: Zolotov et al. . 20141) . The amount of relaxation was shown 
to be sufficient to form classical bulges with high Sersic indices at the center of 
exponential disks ( Ceverino et al. . 20141) . Yet, strong subsequent evolution can oc¬ 
cur and even basic parameters like the bulge:disk mass ratio can largely evolve be- 













































22 


F. Boumaud 


tween high-redshift u nstable phases and present-day galaxies (IMartig et all 12012h . 


Bekki & Cioni ( 2007h highlighted some possible signatures of the past presence of 


giant clumps and contribution to present-day bulges. 

Thus, a consensus on the resulting bulge properties (mass and type) is far from 
being reached, recent efforts having mostly focused on understanding the nature of 
the giant clumps and their own evolution with respect to star formation and feed¬ 
back. Note also that the central coalescence of clumps, while it may induce enough 
relaxation t o produce classical bulges, can significantly e rode a pre-existing dark 
matter cusp ( Elmegreen et al. , 2008 ; Inoue & SaitohL 2011 ). 


3.6 Associated thick disk growth 


Another interesting mechanism associated to the instability of high-redshift disks 
is that pre-existing stars and stars that formed outside the giant clumps or have 
left the giant clumps, are rapidly scattered vertically by the local gravitational po¬ 
tential wells associated to the clumps themselves. The scale-height of the stellar 
disk thus rapidly increases into a very thick stellar system (> 1 — 2kpc). Even if 
the thin disk mass doubles by redshift zero and th e old thick disk shrinks back by 
gravitational response (see Villalobos et al. 2010h . a thick disk will remain with a 
typical scale height of 500-1000pc ; furthermore this instability-induced thick disk 
is decoupled from the younger thin stellar disk formed at lower redshift; this thick 
disk is a separate comp onent in vertical profiles rather than just a low-density tail 


(IBournaud et alil2009h . 


An interesting property of the thick disks formed through this instability mech¬ 
anism is that their gr owth is concomitant to bulge growth, po ssibly accounting for 
chemical similarities ( Chiappini et al. . 20091 : Chiappini . 2009). Another noticeable 
property is the relatively constant thickness of the thick disk throughout the radial 
profile. While this fails to account for the outer flaring observed for thick disks, 
which is probably better explained by minor merg ers and distant tidal interactions 


(I Villalobos & HelmiL 120081: iDi Matteo et al.Ll201 Ih . it does successfully account for 


the presence of a th ick disk in the innermost reg ions, around central disks and 
bulges, as observed (IDalcanton & Bernsteinl 120021) . This latter property could not 
be explained by minor mergers and tidal interactions, w hich stir and thicken prefer¬ 
entially the low-densi t y outer regions of the stellar disk (IVillalobos & HelmiL 12008 


Bournaud et all 2009; Martig et ah . 2012; Di Matteo et al. . 2011 ). At the same time, 
clump instability cannot account for all observe thick disk properties and a contri- 
buti on of other processes s uch as minor interactions and mergers is likely required, 
too ( Inoue & Saitoh . 2014 ). 

Therefore, while interactions and mergers appear needed to explain the outer 
structure of thick disks, clumpy disk instabilities at high redshift appear required 
to explain their inner one. An interesting property of thick disks is that the frac¬ 
tion of the stell ar mass that they gather is la rger in later-type galaxies with small 
bulge fractions (lYoachim & DalcantonLl2006h . This relation might be explained by 
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clump-driven bulge growth and disk thickening, if the early-formed stellar mass 
is distributed between the bulge and the thick disk, with preference for the bulge 
at high total mass and preference for the thick disk at lower total mass. The rela¬ 
tion between bulge properties and thick disk properties, potentially resulting from 
the role of high-redsh ift disk instabilities in bulge growth, was further outlined by 


Comeron et al.l (120141) who argue for concurrent growth of the thick disk and central 


bulgy mass concentrations in the past history of today’s spiral galaxies. 


3.7 Is bulge formation too efficient? Stellar feedback and bulge 
growth regulation 


A key question related to bulge formation or growth by disk instabilities is whether 
this mechanism would over-predict bulge formation. Standard A-CDM galaxy 
formation models already tend to over-produce bulges and spheroids, at the ex¬ 
pense of high angular momentum disks, especially in the baryonic mass range of 
lO*** — 10*^ Mq. At best, some baryonic physics models may resul t in an accept¬ 


able d istribution of stars among bulges and disks (lA^ertz_e^d 


201 ih but generally over-produce the stellar mass dGuo et al.L 


2011: Guedes et al 


201 ll) . while models 


with more realistic stellar masses remain too dominated by bu l ges an d low angular 
momentum components ( Scannapieco et al. . 2011 : Guo et al. . 201 ih . While these 
results are often in tension with obseryations, or can at best be marginally recon¬ 
ciled, they are generally supported by semi-analytic models that do not include disk 
instabilities, or onl y low-redshift secular in stabilities (bars) that grow bulges much 
more slowly (e.g., ISomeryille et all l2008h . or cosmological hydrodynamic simu¬ 
lations that generally do not resolye giant clumps and other features associated to 
disk insta bilities, as this require s high resolution that is too costly to maintain down 


to z = 0 ( Ceyerino et al. . 201 (lh - some hydrodynamic models also employ ther¬ 


mal models and/or feedback sche mes that suppress strong di sk instabilities, as these 
were not obseryed until recently ( Sorneryille & Da\^ 2014 ). Note also that cosmo¬ 
logical simulations tend to oyerproduce stars at earl y epochs (z > 3) and pres erye 
too low gas fractions down to redshifts z = 1 — 3 (iDekel & Mandelkeiil2014l and 
references therein) which can also damp the disk instability process in these simu¬ 
lations. 

The new mechanism of bulge formation by disk instability thus comes on top of a 
cosmological model which, depending on the assumed (and still uncertain) baryonic 
physics, already produces enough stellar mass in bulges and low angular momentum 
components - if not already too much. This could be an indirect argument against 
clump suryiyal and coalescence into bulges, although the obseryed clump ages are 
consistent with long lifetimes and migration to bulges. Note neyertheless that the 
instability driyen inflow (see aboye) is independent of clump suryiyal so its con¬ 
tribution to the growth of central compact components could not be suppressed in 
the case of short-liyed clumps. The question of whether the proposed high-redshift 
disk instability mechanisms oyerproduce bulges is thus naturally raised. The typi- 
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cal numbers for a galaxy of stellar mass about 5 x IO'^'Mq are, say, five clumps 
of 5 X 10^ Mq each, migrating to the bulge in 400 Myr, with unstable steady state 
maintained for 2Gyr. This means that 25% of the total stellar mass coalesces into 
the bulge through this process, leading to an excessive bulge:disk mass ratio of 1:3. 
Furthermore the diffuse instability-driven inflow may double this estimate, which 
was based based on the migration and coalescence of giant clumps only. The result¬ 
ing bulge fraction could thus be above the acceptable levels for such moderate mass 
galaxies, and even if the thin disk doubles its mass without further growth of the 
bulge between z = 1 — 2 and z = 0, the resulting bulge:disk ratio at z = 0, of at least 
1 :6, still seems in tension with observations - especially if other processes, such as 
minor mergers and some bigger interactions that are unavoidable at some level, also 
grow the central bulge. These simple estimates highlight the potential issue. 


Early models of clump formation and migration clearly over-produced bulge 
masses, with final bulge:disk mass ratios about 1:1 after the violent instability pe¬ 
riod (e.g., Noguchi , 1999t Bournaud et al. . 2007 ). A strong limitation of these mod¬ 
els was the lack of stellar feedback other than weak supernovae feedback. While 
the global star formation rate of the galaxies were somehow regulated to realistic 
values on the Main Sequence, the clumps did not produce gaseous outflows at re¬ 
alistic rates close to their own star formation rates. As a result, not only they keep 
all their initial mass, but they accrete surrounding material and their cross-section 
for capture of new material only increases. While the initial clump mass in these 
models is in agreement with observations (with typical masses of a few 10 ^ Mq for 
the main few clumps, and rarely more extreme cases) during the first 50-100 Myr, 
their masses eventually become excessive with frequent clumps above 10 ® Mq, and 
the resulting bulge masses after clump coalescence are too high as well. 


The realisation that clumps actually have their mass content regulated by feed¬ 
back, with outflows compensating for the sustained gas accretion, helps solve the 
problem in two ways. First, the clump masses are regulated to a value fluctuat¬ 
ing around their initial mass throughout their migration in the disk, reducing the 
mass available for bulge coalescence by a factor of a few. Second, the fact that 
clumps gradually lose their aged stars and re-accrete gas means that they remain 
gas-rich and even gas-dominated throughout their lifetime (Bournaud et al., 2014 - 
without these processes they would accumulate large amounts of stars and become 
star-dominated before reaching the bulge region). More than half of the material 
reaching the bulge region is thus gaseous, and will form a rotating disk component, 
which will gradually turn into a stellar component later on and could be scattered 
into a bulge component subsequently, but also remains sensitive to feedback pro¬ 
cesses and may be ejected from the central regions through stellar or AGN-driven 
outflows. Detailed simulations quantifying the bulge growth for long-lived clumps 
having their mass regulated by supernovae, photo-ionization and radiation pressure 
feedbacks, show that the growth of bulges by clump migration and instability-driven 
inflows remains at fully acceptable levels, at least after 0.5-1.0Gyr of evolution, for 
galaxies of baryonic mass of 10*®^** Mq (see Bournaud et al. 2014 for quantitative 
results). 
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Whether a long-lasting unstable steady state will eventually over-produce bulges, 
or whether feedback processes can prevent too much gas to accumulate in the central 
kpc and turn into stars that will subsequently be scattered into a spheroidal bulge, 
remains an open question. Stronger regulation than just clump outflows appears 
needed if the clumps are actually long-lived and can migrate. A first solution to this 
problem was proposed by Ferret et al. (2013). As clumps remain gas-rich along their 
evolution, their central coalescence conveys large amounts of gas inwards (typically 
a few 10* Mq of gas over ^0.1 kpc^) which provokes a local starburst in the central 
kpc. The resulting feedback expulses large amounts of gas - from the coalescing 
clumps, as well as gas brought inward by the global instability-driven inflows. Gas 
expulsion rates from the central kpc peaking at a few tens of solar masses per year 
are reported by Ferret et al. In addition, this process will affect any pre-existing stel¬ 
lar bulge. The gas mass indeed represents roughly half of the mass in the central kpc 
(dark matter providing only a minor contribution at such scales) and this major com¬ 
ponent rapidly fluctuates by clump inflows and feedback-driven gas expulsion. The 
orbits of stars in the bulge are affected by t he rapid fluctuatio ns of the gravitational 
potential - a process already pointed out by Bois et al. ( 2010h regarding stellar scat¬ 
tering by young star clusters and gas clouds in mergers, but involving much higher 
masses at high redshift. 

This process of bulge self-regulation by the inflow of giant clumps is illustrated 
in a simulation similar to those in Bournaud et al. (2014) in Figure |6] In a galaxy 
of stellar mass 5.3 x IO'^’M© with 43% of gas during the analysed period, a giant 
clump of 3.6 X 10* Mq of gas and 2.3 x 10* Mq of stars coalesces with the central 
bulge, the stellar mass of which is 3.3 x 10^ Mq. The clump brings gas to the central 
kpc, and young stars (from the clump stellar content and from central star forma¬ 
tion in the clump gas). The local starburst consumes about 60% of the clump gas 
within 25Myr, with a star formation rate in the central kpc peaking at 13Myr in 
Figure |6] Note that this is a high surface density of star formation rate in the central 
kpc, but it does not drive the entire ho st galaxy outside of the typical scatter of the 


Main Sequence (ISchreiber et al.Ll2014ll and the entire host galaxy does not turn into 


a starburst. The gas mass in the central kpc then decreases by a larger amount than 
just gas consumption by star formation, under the effect of stellar feedback-driven 
outflows. The response to the sudden mass increase (clump accretion) and decrease 
(gas outflows) affects the stellar content of the bulge, with some aged stars leav¬ 
ing the bulge region toward an extended stellar halo on high-eccentricity orbits. In 
the particular example shown in Figure |7] the bulge gains about 2x10* Mq from 
the clump coalescence process, but rapidly looses 2.3 x 10* Mq of aged stars, thus 
slightly reducing the total bulge mass after this clump coalescence eve nt. Detailed 


statisti cal studies remain to be performed, but the simulations studied in iFerret et al 


( 20141) show that the bulge mass can be regulated to reasonable amounts that do 
not exceed 10-15% of the stellar mass for Milky Way-mass galaxies, owing to the 
three regulation processes listed above - 1: regulation of the clump mass by steady 
outflows and dynamical loss of aged stars, 2: regulation of the gas richness of the 


“ After applying Gaussian time smoothing of FWHM 2 Myr to erase fluctuations related to the 
numerical sampling of star formation. 
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Fig. 7 Bulge evolution during the central coalescence of a giant clump (simulations from Boumaud 
et al., 2014). The panels display mass of gas in the central kpc (radius 500 pc, top), the mass of 
stars younger than 100 Myr plus the current time (i.e. younger than 100 Myr at f = 0 and younger 
than 150 Myr at f=50 Myr, middle), and the mass of stars older than 200 Myr plus the current age 
(bottom, still in the central kpc). A massive clump coalesces with the bulge at t ~ 60 Myr. The 
gas mass increases when the clump comes in, and decreases due to feedback-driven outflows. The 
mass of young stars in the bulge increases, but that of old star decreases, and the bulge mass is 
regulated to an almost constant value (even slightly decreasing in this example, although its Sersic 
index increases - see text for details). This bulge regulation mechanism was proposed by Ferret 
and collaborators (PeiTet et al., 2013, and Ferret, FhD thesis, 2013). 
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clumps by re-accretion of gas from the disk, and 3: regulation of the central bulge 
mass by central starbursts and relaxation during clump coalescence. A detailed ac¬ 
counting in cosmological context however remains required to study the bulge mass 
budget over the long-lasting clumpy unstable state fr om z > 3 t o z 1. 

An alternative solution was recently proposed by CombesI ( 2014 ) who demon¬ 
strated than in MOND dynamics, realistic clumpy disk are still predicted by sim¬ 
ulations of z = 2 galaxies, but the efficiency of bulge formation is lowered as the 
clump migration timescale increases. Gravitational torquing and inflows should still 
be present, but without central relaxation through clump coalescence, they may form 
only a pseudo-bulge. 


4 The associated growth of supermassive hlack holes 


The violent instability of high-redshift galaxies brings large amounts of gas toward 
their central regions. This is achieved through the migration of gas-rich clumps, and 
more generally by a global inflow of gas driven by gravity torques between the dense 
features arising from the instability, and compensating for the turbulent losses. Sim¬ 
ulations of this process (Bournaud et al., 2011) have shown that an inflow of about 1 
solar mass per year persists down to the central few parsecs, as the gas is not entirely 
depleted into star formation. It is then sufficient to have one percent of the mass 
brought to the central pc accreted by the central supermassive black hole (SMBH) 
to grow this SMBH in realistic proportions compared to the usual scaling relations 
for bulges and SMBHs. Various small-scale mechanisms in the central parsec can 
indeed lead 1% of the available inflowing gas mass to be accreted by the SMBH 
(ICombesl 1200 ih . The process was studied in detail in iGabor & BournaudI j2013h 
who has shown that bright Eddington-limited episodes of Active Galactic Nuclei 
(AGN) accretion can be triggered by the disk instability, and could contribute to 
the bulk of the supermassive black hole mass growth at z = 1 — 3 for galaxies of 
stellar mass lO'**^** M©. The process is illustrated in Figure|8] In the broader cos¬ 
mological context, the role of cold gas accretion onto gas-rich galaxies and internal 


instabilities was probed bv iDubois et al. l(l2012l) : lDubois et al.l (I2013h . 


Observationally, there is a general lack of correlation bet ween the occurrence of 
AGN and morphological signatures of major mergers (e.g.. iKocevski et al.ll2012h 
except in the most luminous QSOs found preferentially in major mergers. In fact, 
moderately bright AGN that drive the bulk of S MBH growth are mostl y located in 
normally star-forming. Main Sequence galaxies ( Mullanev et al.L 2012 ). which are 
generally clumpy unstable disks at z = 1—3. Searches for a direct link between disk 
instabilities and AGN are hampered by the high gas column densities, typically a 
few times 10® Mq kpc^^ in the central regions of these galaxies, sufficient to reach 
Compton thickness or at least severely attenuate the X-ray signatures of potential 
AGN, indep endently of their feedin g mechanisms. Using optical line emission to 
probe AGN, iBournaud et al.l ( 2012 ) haye shown that at intermediate redshift (z ~ 
0.7) there are more AGN in clumpy unstable disks, than in regular smooth disks 
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Fig. 8 AMR simulations of a high-redshift gas-rich disk galaxy with a central SMBH, using grad¬ 
ual zooms in the AMR refinement toward the central black hole, with a spatial resolution reaching 
0.02 pc in the innermost regions. The white circle in the last zoom represents the SMB FI position 
and its Bondi radius. The disk instability drives steady gas inflows toward the SMBH. AGN feed¬ 
back triggers hot winds that escape through low-density holes and leave the accreting channels 
almost unaffected. The large-scale star-formation activity also remains unaffected in spite of the 
efficient AGN-driven outflows. Figure courtesy of Jared Gabor. 


of the same mass and size (both clumpy unstable disks and modern spira l types 
co-exist at such intermediate redshifts). At high redshift z > 1, Trump et al. ( 2014ll 
find an AGN frequency as high in clumpy disks as in compact early-type galaxies, 
which are known to be frequent AGN hosts (compared to star-forming spirals at 
low redshift), confirming the efficient feeding of AGN in clumpy disks. Yet a direct 
comparison of AGN feeding in clumpy unstable disks and in more “stable” disks is 
impossible at z > 1, as stable spiral disks are virtually inexistent at z > 1. 

High-redshift disk instability can contribute to SMBHs in two other ways. First, 
the clumpy accretion onto black holes can help increase their mass more rapidly at 
early epochs, which subsequently increases the limiting Eddington rate and makes 
possible for the SMBH to grow its mass more rapidly. This potential solution 
to the problem of ve ry massive and bright AGN a t very high redshifts (z ~ 6, 
Di Matteo et al. 2012 ) was studied in DeGraf et al.l ( 2014 ). Second, the clumps 
could be the formation site of SMBH seeds, if they form intermediate mass black 
holes through runaway stellar collisions, which their estimated star formation rate 
densities make possible. These seeds could be gathered centrally along with clump 
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migration into an SMBH (lElmegreen et al.Ll 2008h. w hich w as potentially supported 


by some observed spectral signatures (IShaniro et al.L 120091) . 


Studies of the response to feedback show that the feeding of AGN by disk in¬ 
stability does not quench star formation, and not even the fuelling of the AGN it¬ 
self. AGN from clumpy disks produce high-velocity outflows that are collimated 
by the density and pressure gradients in the disk, and escape perpe ndicularly from 
the disk plane from the nuclear region ( Gabor & Bournaud . l2014l) . This is agree¬ 
ment with recent obs ervations of high-veloc i ty win ds emerging preferentially from 


the nuclear regions (IForster Schreiber et al.L l2014l) . Outflows from star formation 


are more widespread abov e the entire disk and its star-form ing clumps (observa- 
: Newman et al.l2012 , simulations; Bournaud et al.l20T^ Hopkins et al., 2013). 


tions; 


Hence the AGN feedback does not affect the inflowing, which is the fuel for future 
AGN feeding, and the extended gas disks including its star-forming r egions. This 


holds even once long-range radiative effects are taken into account (IRoos et al 


20141) . The AGN luminosity and accretion rate strongly fluctuate over Myr-long 


timescales, which results from the high heterogeneous, turbulent nature of the in¬ 
flowing gas, rather t han from the regulation by feedback ( Gabor & Bournaud. 2014 ; 


DeGraf et al.L 120141) . 


5 Disk instabilities and early-type galaxy formation 


The instability-driven inflow scales like the circular velocity squared (Section 3.4) 
and is thus much more intense in high-mass galaxies. Clump migration is also faster, 
following the dynamical friction timescale in massive galaxies with high-density 
disks and halos. This raises the question of whether the violent instability of high- 
redshift galaxies can form early-type galaxies (ETGs) at high masses, i.e. entirely 
spheroid-dominated systems rather than just bulges in the center of disk-dominated 
systems. 

Theoretically, these strong inflows can lead to disk contraction in a timescale 
not larger than 1 Gyr, and the instability-driven bulge growth rate could lead to a 
bulge-do minated system, through which the disk is stabilised and star formation is 
quenched ( Dekel & Burkerd l2014l). with the help of stellar spheroids stabilising gas 
disks to quench star formation JlVlartig et al. . 20091) . Recent cosmological simula¬ 
tions have probed these possible mechanisms, where the strong inflow first forms 
wet compact star-formi ng systems, which are subsequently quen ched and turned 


into red compact ETGs (ICeverino et al.Ll2014tlZolotov et al.Ll2014l) . The high Sersic 


indices and dispersion-dominated kinematics are consistent with these being the 
progenitors of modern ETGs. The transitions from the compact star-forming system 
to a quenched one could correspon d to observations o f the so-called “blue nuggets” 


and “red nuggets” at high redshift (IBarro et al.Ll2014l) . 


Observations of giant clumps or clump re mnants in the innermo st regions of 


young E TGs in the Hubble Ultr a Deep Field (lElmegreen et al.L 120051) support this 


scenario. iBournaud et al.l (1201 ll) have also shown that mergers of gas-rich unstable 
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disks lead to compact spheroid formation when the instability in the cold interstellar 
phase is taken into account during the merger. Nevertheless, it remains unknown 
whether these processes can explain the detailed phase space stru cture of modern 


ETGs , including the observed families of fast and slow rotators (lEmsellem et al 


20071) which could also be relatively well explai ned in the cosmo logical context 


without invoking a major role of disk instabilities (iNaab et al.L 120141) . 


6 Comparison to secular disk instabilities at lower redshift 


The violent instability where the entire disk is self-regulated at Q 1 persists until 


about redshift 1, for the galaxy masses that we have studied here (lElmegreen et al 


2007t Genel et HI 2012 : Dekel et al. . 2009bt Ceverino et al.L 2014). This violent 
phase, with irregular clumpy disks, growing spheroids, and relatively frequent merg¬ 
ers, has a morpholo gy poorly correlate d to the final bulge/disk ratio of today’s 
descendent galaxies ( Martig et al.L 2012h . After z ~ 1, galaxies enter their secular 
phase where a stable thin disk grows and slowly eyolyes, with a bulge/disk ratio 
close to the final yalue. This regime differs by haying globally Q > 1, with Q<\ 
only locally, for instance for gas compressed in s piral arms and in w hich small 
molecular clouds form by yarious local instabilities ( Renaud et all 2013h . The eyo- 
lution from the early yiolently unstable phase to the secular stable spiral disks is 
shown for two typical cases of zoom-in simulations in cosmological context in Fig¬ 
ure |9] 

The mild instability of modern spirals differ from the global instability of their 
high-redshift progenitors in yarious ways. Present-day disks globally deyelop only 
non-axisymmetric (m > 1) modes such as spiral arms and bars, and grayitational 
collapse at 2 < 1 can occur only locally in small oyer-densities of gas. The as¬ 
sociated inflows are much slower, with only a few percent of the angular mo- 
mentum transferred outwards per rotation perio d eyen in strongly barred galaxies 
( Combes & Gerin . 1985 : Bournaud et all 20051) . Mass inflows and yertical reso¬ 
nances can_secufarly grow central spheroids. Yet, simulations in cosmological con¬ 


text in lKraljic et al.l (12012h suggest that the contribution of these low-redshift secu¬ 
lar instabilities in the bulge mass budget is also more modest than for high-redshift 
instabilities, although the secular phase last longer and sometimes grows massiye 
peanut-shaped bulges. The absence of yiolent relaxation, unlike the central coales¬ 
cence of giant clumps, is such that the process mostly results in pseudo-bulges rather 
than classical bulges with high Sersic indices (Chapter 6.3). 

The inflow toward AGN and SMBH is also much more modest, but another 
difference here is the common presence of Inner Lindblad Resonances (ILRs) in 
low-redshift spiral galaxies - high-redshift disks generally haye no ILR associated 
to the clump instability. As a consequence, the inflow stops and the gas is stored 
at the ILR radius until a nuclear instability occurs and brings the gas reseryoir 
inwards. The process can be repeated with cyclic AGN feeding, but a large frac¬ 
tion of the inflowing gas can also be depleted through star formation in the mean- 
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Fig. 9 Simulations in cosmological context, zoomed on individual galaxies (from Martig et al., 
2012), displaying the stellar mass surface density (panel size: 20x20 kpc). These simulations show 
the transition from a “violent phase” at z > 1 with violent disk instabilities (V.D.I., top) and giant 
clumps and, more rarely, merger-driven starbursts, to a “secular phase” with bars and spiral arms at 
z < 1. Interestingly, this transition shows a “downsizing” behaviour with stellar mass, i.e. it occurs 
later-on for lower-mass galaxies, which could explain that clumpy disk instabilities persist longer 
for lower-mass galaxies (Elmegreen et al. 2007, Boumaud et al. 2012) and regular barred spiral 
morphologies arise earlier-on for high-mass galaxies (Sheth et al. 2008, Kraljic et al. 2012). 


while (lEmsellem et al.L 120151) . Indeed, observations point out that the correlation 


between galactic bars an d AGN is complicated by a number of factors (see e.g. 
Coelho & Gadottii 2011 and references therein). However, it can b e argued that th e 


correlation between nuclear bars and AGN is more straightforward (ICombesil200lh . 


7 Summary 


High-redshift star-forming galaxies at z ^ 1—3 have irregular optical morphologies 
dominated by a few bright giant clumps, also faintly detectable in the near-infrared. 
There is broad evidence that these giant clumps (of a few 10^^® Mq of gas and stars 
and 500-1000 pc diameter for the biggest ones) form mostly by in-situ gravitational 
instability in gas-rich, turbulent galactic disks. This is largely supported by photom¬ 
etry, kinematics, and stellar population studies. Recently, a first example of direct 
gravitational collapse of a giant clump could be directly probed (in the form of a 
very massive collap sing star-forming blo b with only very limited and very young 
stellar counterparts, Zanella et al.L 2015 ). Only a small fraction of clumps exhibit 
older stellar populations and may form ex-situ, in the form of small satellites of 
gaseous clumps that merge with the disk from the outside. 
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The modern understanding of galaxy formation in the standard cosmological 
framework explains the high gas fractions and resulting disk instability as the out¬ 
come of steady accretion of cosmological gas reservoirs (and some companion 
galaxies) at high mass rates, keeping the gas fraction high, the Toomre stability 
parameter low, and the disk in a globally unstable state. The disk increases its gas 
velocity dispersion (or turbulent speed) to self-regulated its dynamics in a steady 
state about 2 ~ 1. 

Hence the high-redshift progenitors of Milky Way-like spirals differ from mod¬ 
ern disk galaxies, which have only weak non-axisymmetric instabilities (bars and 
spiral arms) and in which gas undergoes gravitational collapse only in limited re¬ 
gions, in the form of transient low-mass molecular clouds. 

The giant clumps and the underlying instability can build a galactic bulge in sev¬ 
eral ways. The first one is the instability-driven inflows, which pumps gravitational 
energy into the interstellar turbulence cascade to compensate for the radiative losses. 
This inflow builds a central mass concentration in the form of a pseudo-bulge - un¬ 
less another process increases the relaxation and turns this central concentration into 
a classical bulge. 

Detailed numerical models of star formation and feedback in a multi-phase ISM, 
and observations of stellar population ages, mostly support that the giant clumps 
can survive against feedback from young massive stars for a few hundreds of Myr, 
unlike nearby molecular clouds. In this case the giant clumps undergo dynamical 
friction from the host galaxy and its dark matter halo and migrate inward in a few 
10^ yr, coalesce with the central bulge, or form a bulge if no bulge is present yet. In 
this case the induced relaxation is generally found to turn the central spheroid into 
a classical bulge with a high Sersic index. 

The detailed properties of bulges built by disk instabilities remain uncertain, and 
highly dependent on the physics of stellar feedback. Some studies find that it could 
be possible to form only a pseudo-bulge, but others And that in the most massive 
galaxies the whole system may turn into a classical spheroid through the violent 
disk instability, consistent with properties of early-type galaxies. 

An interesting property of clump migration and central coalescence is that the 
induced relaxation can affect the stellar orbits of a pre-existing bulge, and cause 
dynamical evaporation from the central bulge toward a very extended, low-density 
faint stellar halo. In this case the process gets self-regulated and the bulge mass is not 
grown beyond 10-15% of the total stellar mass for a Milky Way-mass high-redshift 
galaxy. 

The violent instability of high-redshift disk galaxies presents other interesting 
properties that can form other sub-galactic components concurrently with bulges. 
The instability-driven inflow can typically provide one solar mass per year toward 
the central parsec, which may be sufficient to dominate the feeding of central su- 
permassive black holes in moderate mass galaxies, and this is potentially supported 
by observations of active galactic nuclei in Main Sequence galaxies. Along with 
bulges and central black holes, the violent instability of high redshift galaxies can 
also grow the old thick stellar disks, which are ubiquitous around present-day spiral 
galaxies. 
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